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K* ' Abstract. We investigate the effect of resonances in photoionization and 

(N ■ 
t^ . 

^^ , and recombination rates derived from high-resolution cross sections for oxygen ions 

\^ ' are compared with earlier works with less resolution and accuracy, such as in the 

>-p . widely used Opacity Project data. We find significant differences in photoionziation 

^— N I rates for O II metastable states, averaged over Planck functions corresponding to 

"^ ' ionizing radiation fields, with respect to the intrinsic accuracy of the calculations 

^ , and improved resolution. Furthermore, for highly charged ions other physical effects 

I ' are also important. Recombination rate coefficients, averaged over a Maxwellian 

; I ■ distribution, are extremely sensitive to the position and resolution of near-threshold 

j/3 I resonances, and radiation damping, in (e + O VII) < — > O VI + hv. Surprisingly 

^ ' however, the effect on the monochromatic and the mean Rosseland and Planck bound- 



free opacities is relatively small, but may be potentially significant. 



1. Introduction 

A precise treatment of resonances has always been of prime concern in coupled 
channel calculations for electron impact excitation, photoionization, and recombination. 
Sophisticated methods have been employed to resolve and/or fit resonances in order 
that cross sections and averaged rates, the quantities of interest in applications, may 
be computed accurately. One of the main reasons that the powerful R-matrix method 
(Burke et al 1971, Hummer et al 1993, Burke and Berrington 1993) has been very 
successful is because it enables efficient calculations at a large number of energies to 
affect such a resolution. Following a single diagonalization of the R-matrix for each 
total electron- ion symmetry LSvr or Jtt, cross sections may be obtained at any number 
of energies. The R-matrix method, and its relativistic variant the Breit-Pauli R-matrix 
method (BPRM, Berrington et al , 1995), have been employed in large-scale calculations 
in the Opacity Project (OP, Seaton et al 1994, The Opacity Project Team 1995, 1996), 
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and the Iron Project (IP, Hummer et al 1993) for radiative and collisional cross sections 
for most astrophysically abundant elements, nearly 200 atoms and ions. The OP/IP 
work involves photoionization and electron scattering cross sections, as well as transition 
probabilities. The R-matrix method (and BPRM) has been further extended to a 
unified treatment of electron-ion recombination including radiative and di-electronic 
recombination (RR and DR) in an ab initio manner, which is also self-consistent with 
R-matrix photoionization cross sections (Nahar and Pradhan 1992, 1994a; Zhang and 
Pradhan 1997, Zhang et al 1999, Nahar and Pradhan 2003a). A crucial aspect of the 
unified method for (e + ion) recombination is the high resolution imperative for the 
derivation of accurate rate coefficients. However, some physical effects are also very 
important in determining the accuracy of computed parameters as we demonstrate in 
this paper. The present work also addresses some of the issues and results described by 
Ramirez and Bautista (2002). 

2. Results 

The effect of resolution and accuracy of resonances depends on physical processes. 
The height of photoionization resonances is theoretically unbounded, unlike electron- 
ion excitation cross sections where it is limited by the unitarity constraint on the 
scattering matrix. Therefore resonances in photoionization cross sections may approach 
unphysical heights if all relevant physical effects are not carefully considered. A 
complete resolution of resonances in photoionization cross sections, particularly for 
higher n and i, is difficult to achieve. Neither can the calculated resonances be 
individually verified by experiments, since the measured cross sections are averages over 
the beam width and shape. Nevertheless, R-matrix calculations for photoionization 
cross sections have been experimentally verified, some computed in LS coupling prior to 
experiments. For example, we note the extremely detailed comparison between R-matrix 
cross sections and recent experiments on synchrotron based light sources reported by 
the Reno/Berkeley group for O II (Covington et al, 2001), the Aarhus group for C II 
(Kjeldsen et al. 1999, Nahar 2002), and the Paris group for O ions (Champeaux et al 
2003, Nahar 2004). Near-threshold resonances need to be carefully resolved, often with 
relativistic fine structure. In the relativistic Fe XVII photorecombination calculations 
including 359 bound levels with n < 10, Zhang et al (2001) explicitly delineate 2,985 
resonances in a very small near-threshold region - 5eV. 

2.1. Photionization cross sections and rates 

In Fig. 1 we compare high-resolution R-matrix photoionization cross sections 
for the ground and two metastable states of O II, 2s'^2p^{*S°,'^ D°,'^ P°), with 
the OP data from Topbase (Cunto et al 1993, http://vizier.u-strasbg.fr/OP.html, 
http://heasarc.gsfc.nasa.gov/topbase/home.html). The OP cross sections were the most 
accurate attempt at the time of computations to resolve resonances. In general, instead 
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of a constant energy mesh, a quantum-defect mesh was employed to enable similar 
resolution of resonances along a Rydberg series as they get narrower with increasing 
n. The quantum-defect mesh ensures a pre-determined number of points (usually 100) 
to resolve each range of effective quantum number (z/, u + 1). Although this does not 
guarantee complete resolution, and the OP cross sections may still be under-resolved 
for certain applications such as monochromatic spectral modeling, they are considerably 
better than those with constant energy mesh. 

The new results reported in this work are obtained as in the recent calculations for 
a number of oxygen ions, O II - O V, using a more extensive eigenfunction expansion 
for the target ion than in the OP work; they are in good detailed agreement with 
recent experimental measurements (Nahar 2003). Extensive resonances structures are 
delineated, especially in the near-threshold region. As Fig. 1 shows, resonances are 
more completely resolved in the present calculations, including those along a Rydberg 
series converging on to excited thresholds. Moreover, one important point is that the 
resonances positions are at somewhat different energies compared to the corresponding 
OP data. 

Photoionization rates Tpj are computed using both the present and the OP cross 
sections, assuming ionizing radiaton fields from a blackbody source {By{T)) at radiation 
temperatures T represented by a Planck function. 

J hu 

B^{T) = {2huyc^)[exp{hu/kBT) - 1]~^ (2) 

where W is a dilution factor (set to 1 for the present study) and apj is the photoionisation 
cross section. 

We can see in Fig.l, the ground state cross sections agree well. The position of 
the resonances are identical. Their heights differ but, as shown below, it should not 
affect the rates significantly because the area under these peaks is not significant since 
the resonances are narrow. For the two metastable states a shift in the position of the 
resonances can be seen. The difference in resolution also modifies the computed shape of 
the resonances, especialy for the first ones. These two factors, resolution and accuracy of 
the position of the resonances, results in significant deviations in photoionisation rates. 

In Fig. 2 we present the percentage difference in the rates, obtained using OP cross 
sections and the present data, for the three lower states of O II at radiation temperatures 
ranging from 10'^ — 10^ K. While there is little differences for the ground state ^5*°, 
because the cross sections agree very well in the position of the resonances as well as 
in the resolution between the two set of data, the differences for the two metastable 
states are quite signficant, up to 30% for ^P°, and up to 55% for the ^D". In order 
to ascertain whether the differences are due to resolution, or accuracy as reflected in 
the different positions of resonances, another set of rates are obtained on shifting the 
resonance positions by a slight amount to match in the two sets of data. This reduces 
the differences by 10%, with new maximum differences of 20% and 45%. Resolution 
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also plays a role in the discrepancy. For the metastable states (^P°, ^-D°), the first 
resonance is much broader in the OP data. The areas under this first peak are 1.8 
and 1.9 times higher respectively for the OP cross sections compared to the present 
set. This affects the low temperature regime more severely since the first resonance is 
crucial in convolving with the Planck function, although less so than with a Maxwellian, 
considered next. 

2.2. Recombination rates 

A similar situation exists for the highly charged ion Li-like O VI. In previous calculations 
it was shown that the large KLL resonances were severly under-resolved in previous 
works, and a very high resolution is required in order to obtain converged values for the 
resonance oscillator strengths. The peak values of resonances rise more that 4 orders 
of magnitude above the background cross sections. Moreover, since the core ion is 
He-like, with a large radiative decay rate for the dominant core ls2j9(^P°) -^ 15^(^5*0) 
transition, radiation damping is important (Zhang et al 1999). This is particularly 
so in the calculation of unified photo-recombination rate coefficients for (e + O VII) 
< — > O VI + hz/ (Nahar and Pradhan 2003b). The unified scheme to obtain total 
recombination rates, including RR and DR, entails: (A) highly resolved calculations 
of photoionization cross sections of a large number of bound levels, typically several 
hundred levels with effective quantum number v < 10, and (B) DR calculations using 
an adaptation of the Bell and Seaton theory (1985, Nahar and Pradhan 1994a) for 
10 < z/ < oo. In addition, for any given ion self-consistency is assured in an ab initio 
manner between total photoionization and recombination since the same wavefunction 
expansion is employed for both processes (see Nahar 2003, and references therein). 
Another advantage of the method is that level-specific cross sections and rates for all 
levels with z/ < 10 are also obtained. 

Total and level-specific rate coefficients for recombination to O VI and O VII 
have been reported by Nahar and Pradhan (2003b). In a recent work, Ramirez and 
Bautista (2002) calculated photoionization cross sections for the ground state of highly 
charged ions such as O VI and Fe XXIV using a method to locate resonances for 
high resolution (Quigley and Berrington 1996), but neglecting effects such as fine 
structure and radiation damping. Therefore their photo-recombination rate coefficients 
are in error by large amounts, as seen in Fig. 3; there is over an order of magnitude 
discrepancy relative to Nahar and Pradhan (2003b). Thus in this case it is seen that both 
resolution and accuracy, as reflected in the neglected physical effects, are important. 
In their work Ramirez and Bautista (2002) discuss the earlier photoionization and 
(e + ion) recombination calculations for Fe XVII by Pradhan et al (2001); however 
their uncertainly estimates of "25-50%" for "integrated resonance contributions" are 
not derived from direct computations, and are much higher than those computed by 
Zhang et al (2001). 
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2.3. Opacity 

Finally, we investigate these effects on a larger scale, in the total opacity of O II. Fig. 
4 shows the monochromatic bound-free (bf) opacity due to O II (in atomic units) at 
a representative stellar temperature of 2 x 10^ K and electron density of 10^^ cm^^, 
computed as described in Seaton et al (1994). 

«:^^(z/)=^iV,af(z/)[l-e-^] (3) 

i 

The top panel shows results using the present high-resolution photoionization data 
for 329 bound states of O II, compared with results from OP data. At this temperature 
and density, 197 bound states of O II contribute to the bf-opacity in the present 
calculations, and 169 in the OP data. The difference in number of bound states is 
inconsequential; the additional states in the present work are mostly highly excited ones. 
The monochromatic opacity is sampled at 10^ points in both calculations. As Fig. 4 
shows, there is overall very good agreement between the two sets, despite significantly 
better resolution and accuracy of the present calculations. This is expected since strong 
features dominate the opacity spectrum, and are found in both sets of calculations. The 
present Rosseland and Planck mean bf-opacities 

1 r 1 



bf - I bfFR^^)^^ (4) 

R '' l^v 

K^J = f kI^ Fp{u)du (5) 



where 



Fni'^) = —Air^r- 73^ (6) 





hu 


ksT 




hu 



- kgT 



Fp{u) = —,{Tr^Y — (7) 

are 0.159 and 2.12 respectively, compared to 0.146 and 2.42 from the OP data, 
i.e. the Rosseland means are 8% higher, and the Planck means are 14% lower than OP 
(electron scattering opacity is included to provide a background, and is the same in both 
calculations). In the present opacity calculations the resonances are better resolved, and 
in many cases narrower compared to the OP data, which results in a smaller Planck 
mean; conversely, since the minima in resonances are better sampled, the Rosseland 
mean is slightly higher. However, the total Rosseland mean opacities differ only by 
2.6%, if we include the same (much larger) bound-bound contribution from lines in the 
two calculations — 0.345 from the present vs. 0.336 from the OP data. 
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3. Conclusion 

In conclusion we note that: 

1. For photoionization and NLTE modeling applications, it may be necessary to 
compute photoionization cross sections with higher accuracy and resolution than in the 
Opacity Project data. Accuracy of resonance positions is crucial in the determination 
of (e + ion) and photoionization rates. A number of self-consistent computations are 
in progress as part of a program to calculate photoionization and unified (e + ion) 
reocmbination (RR plus DR) cross sections and rates (Nahar and Pradhan 2003a). 

2. Although the monochromatic and mean opacities are not measurably different 
using the high-resolution data for oxygen ions, as opposed to the Opacity Project data, 
these quantities may differ considerably for more complex atomic systems such as the 
low ionization stages of iron. In a previous work (Nahar and Pradhan 1994b), we showed 
that the Rosseland mean opacities using new improved data for Fe II differed by up to 
50% from OP values. There is increasing evidence that there is "missing UV opacity", 
possibly due to Fe I, that is critical to the determination of abundances in the Sun and 
other objects (e.g. Asplund 2003). New atomic and opacities calculations for iron ions 
are needed with higher accuracy and resolution. 

PACS number: 34.80.Kw 



Resolution and accuracy of resonances in R-matrix cross sections 7 

Submitted to: J. Phys. B: At. Mol. Opt. Phys. 

Acknowledgments 

This work was partially supported by the US National Science Foundation and NASA. 
The computational work was carried out on the Cray SVl at the Ohio Supercomputer 
Center in Columbus, Ohio. 

References 

Asplund 2003 astro-ph 0312291 {Astron. Astrophys. in press) 

Bell R H and Seaton M J 1985 Journal Of Physics B 18 1589 

Berrington K A, Eissner W, Norrington P H 1995 Comput.Phys. Commun. 92 290 

Burke P G, Hibbcrt A and Robb W D 1971 Journal Of Physics B 4 153 

Burke P G and Berrington K A, Atomic and Molecular Processes: an R-matrix Approach, lOP 

Publishing, Bristol, 1993 
Champeaux J -P, Bizau J -M, Cubaynes D, Blancard C, Nahar S N, Hitz D, Bruncau J, and Wuillcumier 

F J 2003 Astrophys. J. Supp. Ser. 148 583 
Covington A M, Aguilar A, Covington I R, Gharailbeh M, Shirley C A, Phancuf R A, Alvarez T, 

Cisneros C, Hinojosa H, Bozek J D, Dominguez I, Sant'Ama M M, Schlachter A S, Bcrrah N, 

Nahar S N, McLaughhn B M 2001, Phys. Rev. Lett. 87, 243002-1 
Cunto W., Mendoza C, Ochsenbein F., Zeippen C.J. 1993, Astron. Astrophys. 275 L5 
Hummer D G, Berrington K A, Eissner W, Pradhan A K, Saraph H E and TuUy J A 1993 Astr. Ap. 

279 298 
Kjeldsen H, Folkmann F, Hensen J E, Knudsen H, Rasmussen M S, West J B, and Andersen T 1999, 

Astrophys. J. 524, L143 
Nahar S N 2002 Physical Review A 65 5207 
Nahar S N 2004 Physical Review A 69 042714 

Nahar S N and Pradhan A K 1992 Physical Review Letters 68 1488 
Nahar S N and Pradhan A K 1994a Physical Review A 49 1816 
Nahar S N and Pradhan A K 1994b Journal Of Physics B 27, 429 
Nahar S N and Pradhan A K 2003a Self-Consistent R-matrix Approach To Photoionization And Unified 

Electron-Ion Recombination, in Radiation Processes In Physics and Chemistry, Elsevier (in press). 
Nahar S N and Pradhan A K 2003b Astrophys. J. Supp. Ser. 149, 239 
Pradhan A K 2000 Astrophys. J. 545 L165 

Pradhan A K, Nahar S N and Zhang H L 2001 Astrophys. J. 549 L265 
Quigley L and Berrington K A 1996 Journal Of Physics B 29 4529 
Ramirez J M and Bautista M A 2002 Journal Of Physics B 35, 4139 
Seaton M J, Yu Y, Mihalas D and Pradhan A K 1994 Mon. Not. R. astr. Soc. 266 805 
The Opacity Project, Vols 1 & 2, 1995, 1996, compiled by the Opacity Project Team (Institute of 

Physics Publishing, Bristol and Philadelphia) — OP 
Zhang H L and Pradhan A K 1997 Physical Review Letters 78 195 
Zhang H L, Nahar S N and Pradhan A K 1999 Journal Of Physics B 32 1459 
Zhang H L, Nahar S N and Pradhan A K 2001 Physical Review A 64 032719 



Resolution and accuracy of resonances in R-matrix cross sections 



100 



^ 



K= 10 




1000 



J^ 



b 



Up 


1 1 1 1 


' II 




40 


1 


- 


30 


1 


- 


20 


1 


- 


10 


ym 


A^_ 





50 


1 1 1 1 1 1 








40 






- 


30 






- 


ao 






- 


10 


7 




A - 



0.2 0.4 0.6 



1000 



-II 


1 1 1 1 1 1 1 








I 2P° (OP) 


J 




50 








Z 




,.,,,,,_ 


^ 


' 






- 








40 










— 


- 


100 


- 




- 


30 










- 


— 


50 


_ 




J 


30 










~ 




















^ - 




- 







- ■ 


, - 


.^ 






I 


- 


t1 




0,05 


1 


0-2 0,4 


0,6 


p 


UU 


1 


i 


A— JuJ^ 








- 




[ 










1 1 1 1 1 1 1 1 1 






- 



2 



E(Ryd) 



0.5 1 

E(Ryd) 



1.5 



Figure 1. Photo-ionization cross sections of the ground state ls^2s^2p^ {"^S") and two 
mcta-stable states ls^2s^2p^ i^P°j "^D") of O II (from top to bottom). The present 
high-resolution data are on the left and the Opacity Project data from Topbase are on 
the right. The insets show the shift in resonance positions and the effect of resolution 
in the near threshold region. 
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Figure 2. Percentage differences in photoionization rates for the ground state 
ls^2s'^2p^{^S°) and metastable states ls'^2s'^2p^ {^ P° ,'^ D°) of O II, at radiation 
temperatures T^ad- 



Resolution and accuracy of resonances in R-matrix cross sections 



10 



10- 



10 



-12 



10- 



10- 



10- 




10- 



e + OVII->OVI 

Ground state: ls^2s(^S 

Present(Unified) 

Ramirez & Bautista(2002) 



1 1 1 1 I I 



I I 1 1 1 1 1 



10 



100 1000 



10^ 



lOS 



10^ 



10"^ 



T(K) 



10^^ 



109 



Figure 3. Photo-recombination rate coefficient for the ground state of O VI 
(ls^2s(^S'i/2): solid line - using cross sections from Nahar and Pradhan (2003), dashed 
line - Ramirez and Bautista (2002); the differences are due to neglect of radiation 
damping in the latter calculations. 
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Figure 4. Monochromatic bound-free opacity of O II using the present high-resolution 
atomic data for photoionisation cross sections of 329 bound states (upper panel) , and 
using the Opacity Project data for the same states (lower panel). 



